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Abstract 

The production of quarks and antiquarks during a sudden restoration of 

chiral symmetry, as it might occur in very energetic heavy ion colhsions, is 

considered. If gluons are already present they can assist additively to the 

overall production: real gluons are partially decaying during such a chiral 

transition. The total number of produced quarks is calculated and found to 

be quite sizeable. It is speculated that such a phenomenon could give rise to a 

significant contribution to the overall quark pair creation in the preequilibrium 

stage of the heavy ion collision. 
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One primary goal in the upcoming relativistic heavy ion colhsion experiments at the 
Brookhaven Relativistic Heavy Ion Collider (RHIC) and the CERN Large Hadron Collider 
(LHC) is the temporary formation and potential observation of a new form of matter, the 
quark gluon plasma, a deconfined phase of QCD [|1|. It is expected that this phase is also 
accompanied by a restoration of chiral symmetry [Q]. At high enough temperatures the mas- 
sive and confined quasiparticles, the constituent quarks, become bare, undressed quarks with 
current quark masses being much lighter than the constituent quark masses. Such a physical 
picture is realized in model descriptions of low-energy, effective QCD lagrangians like the 
Nambu and Jona-Lasinio (NJL) model @,|]: At large enough coupling chiral symmetry of 
the vacuum is spontaneously broken and exhibits a nonvanishing scalar quark condensate 
(gg) ~ —(250 Me V)^. A similar behaviour is experienced in lattice QCD calculations 0. 
The scalar condensate attributes a scalar selfinteraction to the fermion propagator resulting 
in the presence of constituent quarks at low energies, in line with experience from hadronic 
spectroscopy. For higher temperatures, however, the condensate becomes weaker and van- 
ishes at sufficiently large temperatures (on the order of T ~ 200 MeV) [§,0, which is also 
found in recent QCD calculations [[7|. 

Relativistic heavy ion collisions offer the only possible way to gain more insight into the 
phase structure of nuclear matter at such high temperatures (and densities). However, from 
the onset the reaction during any such a collision proceeds far off equilibrium, so that a 
thermal quark gluon plasma may be reached only after some finite time. According to the 
results of the parton cascade model ||^ the quarks will reach, if at all, a chemical saturation 
only eventually in a central relativistic heavy ion collision at collider energies whereas the 
gluons saturate on a significantly smaller time scale (in about half a fm/c) P,p!0|. This 'hot 



glue' scenario was proposed by Shuryak [Q. The reason for this particular behaviour lies in 
the fact that the quarks are only produced in second order by gluon annihilation and the rate 
for this process is small compared to the corresponding one for gluon production. The decay 
of a single gluon into a quark and antiquark is, of course, kinematically not allowed, as long 
as all partons are treated as effective on-shell particles. In this work we want to propose that 



the occurrence of a chiral phase transition can have some profound influence on the quark 
pair production, if the restoration will happen on a very short timescale. The idea is that 
if indeed a deconfined state (equilibrated or not) will be formed in the reaction of a heavy 
ion collision, the 'undressing' of the quarks should already have occurred. Being probably 
still far from equilibrium this chiral restoration may thus happen rather spontaneously. We 
assume that the scalar quark condensate a ~ (gg) would melt on a timescale smaller than 
any scale set by the dynamics. The decay of the vacuum induces a spontaneous creation of 
quark and antiquark pairs. During the chiral restoration transition the quarks are no longer 
on-shell particles, their spectral function contains a wide off-shell spectrum, thus permitting 
also the decay of a single gluon into a quark-antiquark pair. The decay of both the vacuum 
and a substantial fraction of energetic gluons would accelerate the chemical saturation of 
the quark phase space. 

To describe the effect on quark pair production we construct a simple model for fermionic 
matter. Denoting the time of chiral restoration at t = 0, the light quark masses then are 
functions of time 

rUc ~ 350 MeV fort < 
m(t) = . (1) 

I r/ifc ~ 10 MeV fort > 
We assume that no further residual interactions act among the quarks. While ansatz (|I]) 
may appear quite artificial, it captures the essential spirit of what one would naively call 'a 
transition to the phase with restored chiral symmetry'. The constituent vacuum becomes 
unstable so that during the restoration of the masses occupied states of the negative Dirac 
continuum are partially transferred into the positive continuum. The amplitude for this is 
obtained in a straightforward way within the In/Out formalism. A particle with momentum 
p sitting in the negative Dirac sea is described by the incoming (constituent) wavefunction 
t/'I'^"'' = if^^^ = y9p|(x, t) (the arrow '|' denotes a Dirac-state in the negative continuum, 
while 't 'denotes one in the positive continuum). After the change of the vacuum structure 
this wavefunction goes over into — > v^p^'' = ctpV^pj — Ppfp^j where states with the label 
(b) denote the outgoing (bare) wavefunctions. The number of produced quarks is obtained by 
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the projection squared on the appropriate outgoing states defining the asymptotic particles 



The number of produced quarks is depicted in fig. 1. It is found that the spectrum resembles 
the distribution of a saturated hot quark-antiquark plasma with a temperature between 200 
and 300 MeV, although the distribution does not drop off exponentially, but as a power law 
for large momenta, Np ^ {rric — m;,)^/(4p^). The total number of produced quarks is hence 
linearly divergent and should be regulated by a typical momentum cutoff A^; as in the NJL 
model 01 which corresponds to a separation between the low and high momentum degrees 
of freedom. For Ac — 1 GeV the total produced quark density is ~ 1.5 — 2fm~^ which 
again is similar to the quark densities in a plasma for temperatures between 200 and 300 
MeV. In addition we have also shown in fig. 1 the situation for a bare quark mass of 150 
MeV {rric = 550 MeV), a typical value for the strange quark. Although the production is 
suppressed for small momenta, the total integrated quark density is only minorly affected 
{us ~ 0.65 - 0.8 fm-3). 

The mechanism for this multiparticle production can be understood as follows: In the 
current quark picture, what happens is that the 'hadronic' vacuum state | In) corresponds 
to a complex coherent superposition of current quark states which rapidly go out of phase. 
The time At needed for the vacuum to restore chiral symmetry is expected to be propor- 
tional to the inverse of the typical nonpertubative QCD energy scale or the chiral symmetry 
breaking scale AqJ^,^ ~ A,^^ ~ 0.2 fm/c. A similar estimate is based on the kinetic equili- 
bration time of the gluons at collider energies predicted to be of the order 0.3 fm/c P,p!0|. 
A small but finite duration At suppresses the pair production at sufficiently high momenta, 
because the wavefunction of the quark resolves the smoothness of the transition. To quantify 
this we substitute for the time dependent mass (|^) 



|TT|] and takes here the simple form 



= p;Pp{OIn I SI;^d%''^ I Jn) = I /?p 




(2) 
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*m 'rnc + rrib rric-mb -2\t\/r\ /o\ 

m {t) = sgn(t) (1 - e I 1/ J (3) 



and propose a nontrivial ansatz for the wavefunction 



( ^'«(t)cos(/?(t)f/ ^ 



j-i(e(t)-p-x) ^ 



where a, (y9 and e are real function in time. Inserting the ansatz into the time dependent 
Dirac equation yields the following coupled set of differential equations 

e = -pcos[2a) 



2 sin ip cos ip 

ha = -m*{t) + pcos(2a) f ) (5) 

y 2 sin (y9 cos cos J 

hip = psin(2Q;) . 

With the appropriate initial condition for the state in the negative continuum these equations 
can be integrated numerically and then projected on the outgoing particle state. In fig. 2 the 
spectrum of particles produced in a 'smooth' transition are depicted for various choices of r. 
The total duration of the transition is approximately At ^ 1.5 r. As expected, the particle 
number at larger momenta depends sensitively on the choice of r. However, for r < 0.2 fm/c 
(and thus At < 0.3 fm/c) only the high momentum yield is affected. We conclude that if for 
dynamical reasons, as estimated above, the transition from the broken to the unbroken chiral 
phase occurs rapidly enough, the change of the underlying vacuum structure is accompanied 
by a substantial nonperturbative production of quarks and antiquarks. 

The additional possibility of gluon decay suggests an interesting field theoretical problem 
and can be phrased as follows: What is the number of quarks (and antiquarks) produced in a 
first order decay of already present gluons propagating in an external time dependent scalar 
background field? We assume now that gluons interact perturbatively with the quarks, 
although they are in part responsible for the dynamic quark mass. The formulation we have 
carried out by using real-time Green functions [0,0 is appropriate for non-equilibrium 
studies. The number of produced quarks is contained in the '< '-component of the complete 
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one-particle Green function and is given by the projection on the outgoing state in the 
distant future 

= *,=feoo/ ((-0^?f(l)G<(l,2)(7o^J'j(2))) . (6) 

For our purpose we have to specify the fermionic initial conditions. For this the zeroth 
order Green function Go is determined by the full set of wavefunctions defined in respect 
to the incoming vacuum state and evolving nonperturbatively in time according to the time 
dependent scalar background Hamiltonian i.e. 

G<(l,2) = +.5:^W(l)^a(2) 
pi 

G>il,2) = . (7) 

pi 

A perturbative expansion is diagrammatically straightforward. The first order decay is 
contained in the lowest order self-energy insertion S for the quarks, where the average 
distribution (rik) of the real gluons enters explicitly. The number of produced particles can 
now be calculated [Q. By using some general properties of the self energy and the Green 
functions, the change in the quark occupation number induced by the self-energy insertion 
reads 

[S>(3,4)G<(4,2) - S<(3,4)G>(4,2)] (7o<^S(2))} . (8) 

The expression (H) for the amount of produced (or scattered) particles resembles the typi- 
cal form used in kinetic transport theories [|1^,|T^; in particular the two terms in the inner 
bracket suggest a 'Gain' and 'Loss' contribution. Two reasons, however, make further anal- 
ysis cumbersome. The first, purely technical one lies in the fact that the Green functions 
depend on both time arguments separately due to the time dependent external field. The 
second one stems from the occurrence of a number of infrared singularities. Besides the 
pair production the above expression also incorporates the absorption or emission of gluons 



for the quarks produced in lowest order in the vacuum decay (0). A separation of the pair 
production from these processes is not possible because of quantum mechanical interference. 



In the end one has to show that these singularities cancel each other, as they indeed do 
Finally, (|]) can be split into three contributions [0, 

where the 'Gain' and the 'Loss' term both are not simply the squared of some particular 
amplitude. However, as a numerical investigation shows, the first is mostly positive while 
the second is negative. The interpretation stems from the fact that for example the 'Gain' 
part contains quarks scattered into the momentum state p either by emission or absorption 
of a gluon as well as the decay of a gluon into a quark and antiquark. The exotic process of 
the simultaneous creation of a quark-antiquark pair and a. gluon also contributes. The 'Loss' 
term contains the opposite processes. The cancellation of the infrared singularities among 
both terms is easily understood within this interpretation. The designation of a 'masshift' 
part for the last contribution is based on the known fact that in an equilibrated plasma the 
quarks become dressed and dynamically massive by the interaction with the thermal gluons 
I . Such a masshift should enter as a correction to the zeroth order vacuum decay (Ef) 



An ultraviolet divergence arising from interactions with the perturbative gluonic vacuum 
is interpreted as mass renormalization. The vacuum terms can explicitly be separated and 
should be regarded as the first order vacuum correction to the spontaneous decay. However, 
we are mainly interested in the overall effect due to the presence of real gluons. In the 
following numerical study these vacuum contributions are therefore discarded. 

The distribution of gluons, at the time where the restoration of chiral symmetry takes 
place, is probably far from local thermal or kinetic equilibrium. In particular high energetic 
gluons parallel to the beam axis may be present. However, for our present purpose it is 
convenient and sufficient to parametrize the incoming gluons by a thermal distribution with 
a large temperature (in the range of 700 to 1000 MeV) which includes such high momenta. 
The results for a chosen temperature of Tg = 700 MeV are shown in fig. 3. The QCD 



coupling constant as = g'^/{4:Ti) is taken to be 0.3. For low momenta the quark distribution 
is large and positive and exceeds unity by three orders of magnitude. It then quickly turns 
and remains negative up to momenta of 400 MeV. At larger momenta the distribution 
resembles a not fully saturated thermal quark distribution of nearly the same temperature 
as the gluons. The distribution, however, drops faster at still higher momenta which are 
not shown here. The low momentum behaviour is obviously unphysical and shows that 
pertubation theory is not applicable here. If we integrate the calculated distribution over 
momentum space, it turns out that the behaviour at low momentum is irrelevant. The major 
amount of particles is created at larger momenta. An integration up to still only moderate 
momenta p ~ 1200 MeV explicitly demonstrates that about 12x0.36=4.3 quarks per fm^ 
are produced by the initiated decay of gluons (integrating up to momenta ^ 2500 MeV the 
number changes to 12x 0.84=10.1 quarks per fm^)! Further numerical investigation shows 
that the importance of this effect increases with a higher temperature employed for the 
gluons. It turns out that the nearly coUinear gluons with high momenta are responsible for 
the decay This suggests that such a phenomenon should be of particular importance 
in the early stages of the gluon evolution when the chiral symmetry restoration takes place. 
The effect, however, is small if we take a bare mass m?, of 150 MeV. The distribution is 
shifted towards higher momenta. Integrating up to momenta ~ 2500 MeV the number of 
produced quarks is evaluated as 0.34 fm~^. Consequently, compared to the hght quarks, the 
number of strange quarks produced by gluon decay would be strongly suppressed. 

To summarize, the occurrence of a rapid chiral restoration of the vacuum can lead to 
some particular new and interesting phenomena in very energetic relativistic heavy ion 
collisions. In this letter we investigated the pair production of quarks within a simple 
model illustrating such a dynamical transition. The presence of real gluons assist the overall 
production: During the change of the vacuum structure the direct decay of a gluon into 
a quark and antiquark pair is kinematically allowed. According to our finding indeed a 
significant amount of quarks are produced during this process. Such a phenomenon could 
give rise to a nonperturbative dynamical creation of quarks and antiquarks (as well as 
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entropy) in the very early stages of the heavy ion reaction. It is tempting to speculate 
that the total number of produced quarks during such a transition and during the further 
evolution suffices so that not only the gluonic degrees of freedom are saturated, but also the 
fermionic. The 'hot glue' scenario could turn out to be a hot 'quark gluon plasma' after all. 
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FIGURES 

FIG. 1. The spontaneous production of quarks is shown for rric = 350 MeV, rrife 2± MeV (a) 
and for rric = 550 MeV, = 150 MeV (b). To guide the eye a thermal distribution of massless 
quarks is also drawn (for a temperature of T = 200 MeV (c) and T = 300 MeV (d)). 

FIG. 2. The production of quarks is shown for a 'smooth' transition in time of the vacuum 
(compare text). The parameters given in the figure are to be read from the left to the right. 
rric = 350 MeV and is taken to be nearly zero. 

FIG. 3. The first order corrections on the particle distribution due to the presence of real 
gluons are depicted. The gluons are taken to be distributed with a temperature of 700 MeV. (a) 
corresponds to the parameter set rric = 350 MeV, mi, = 10 MeV, (b) to the same as in fig. 1. 
For a comparison a thermal distribution of massless quarks with a temperature of 700 MeV is also 
drawn (c). 
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